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We discuss properties of the strongly correlated two-dimensional electron gas in Si MOSFETs 
at low concentrations assuming that the electron liquid is close to crystallization. Analogy with the 
theory of 3 He is emphasized. 



I. INTRODUCTION 

Recent experiments fl|-|l0f in high-mobility two- 
dimensional electron systems indicate that in the pres- 
ence of disorder they cannot be described by the conven- 
tional single particle localization theory 111]]. A number 
of unusual features have been reported taking place in the 
dependence of kinetic coefficients on the magnetic field 
H and temperature T . They are observed at low T and 
relatively small electron concentration n. This regime 
corresponds to r s = E p /Ep » 1, where E p — e 2 pp/eh 
is the potential energy per one electron, Ep = p F /2m is 
the Fermi energy, e is the dielectric constant, m is the 
bare electron mass and pf is the Fermi momentum. In 
this article we give a qualitative interpretation of the ex- 
perimental results assuming that the electron liquid is 
close to Wigner crystallization. Our analysis is based on 
an analogy with the theory of 3 He. We would like to 
explain the following experimentally observed facts. 

A. The electron system exhibits a "transition" as a 
function of n from a metallic phase, where the resistiv- 
ity p of the system saturates at low temperatures, to 
an insulating phase, where the resistivity increases as T 
decreases. The value of the critical concentration n c de- 
pends on the amount of disorder in the sample and cor- 
responds to r s — r c s 3> 1. On the metallic side of the 
transition at T = and n w n c the value of the resis- 
tance is of order %/e 2 . 

B. At T = and electron concentration sufficiently 
close to critical, magnetic field parallel to the film Hu sup- 
presses the metallic phase and drives the system toward 
the insulating phase |^,||,pT). Thus the critical metal- 
insulator concentration n^Hu ) increases with Hu . 

In the metallic phase (n > n c {Hu =0)) and at small 
T the system exhibits a big positive magnetoresistance 
as a function of H\\ . This magnetoresistance saturates at 
#U > fff(n) and p(Hf)/p(0) » I i§. 

C. In the dielectric phase the system exhibits a big 
positive magnetoresistance as a function of Hu , which 
saturates at H\\ > H?>. However, at given H\\ > H» 
the conductance of the system exhibits a big negative 
magnetoresistance as a function of the component of the 
magnetic field H± perpendicular to the film |Q. 

D. In the metallic phase at Hu = and T < Ep the 
resistance p(T) significantly increases with temperature. 
The characteristic value of d\np/dT > Ep 1 at small T 



is large and depends on the value of n — n c . 

E. If at H\\ > i?| the system is still in the metallic 
phase (n > n c (Hn)), the T dependence of the resistance 
is much smaller than in the = case (dlnp(iJ|| > 
Hp/dT^E' 1 ) ffl. 

F. In the metallic regime in magnetic fields of order Hfi 
or less the system exhibits a large anisotropy: The mag- 
netoresistance in the parallel field is much larger than the 
magnetoresistance in the perpendicular field . 



II. A COMPARISON OF LOW TEMPERATURE 
PROPERTIES OF 3 HE AND THE ELECTRON 
LIQUIDS. 

We believe that the physics of the strongly correlated 
electrons in Si MOSFETs at r s ^> 1 is quite similar to 
the low temperature physics of 3 He, despite the differ- 
ence in interactions. Indeed, in the absence of disor- 
der, both systems exhibit a quantum zero-temperature 
liquid-crystal transition. In both cases the transition is 
believed to be of the first order. The liquid phase near 
the transition in both systems is characterized by a large 
ratio of the potential to kinetic energy. In the 3 He liquid 
the ratio between the Debye energy Qpi and the Fermi 
energy reaches 80 near the melting pressure [^5|. On 
the other hand, it is known, from numerical simulations 
p2| , that in the two-dimensional case the Fermi liquid- 
Wigner crystal transition takes place at a relatively large 
r s = r c s ~ 38. Therefore, in both systems there are wide 
regions of concentrations where the systems are in the liq- 
uid phases, though the potential energy is already much 
larger than the kinetic one. 

The only qualitative difference between these systems 
is that 3 He crystallizes upon increasing the atomic con- 
centration, while the electron system crystallizes when 
the concentration is lowered. 

A model which describes properties of liquid 3 He 
near the crystallization pressure has been introduced 
in |l5| , |r7| , p^| . It is based on the existence of the strong 
inequality hflp, 3> E F ^ He ^ which indicates that the 3 He 
liquid can be regarded as a system which is being close to 
solid |l^,|T^,^8| . Since the system has strong short range 
crystalline order correlations, on a fast time scale, Ti/ilp, 
the atoms vibrate near certain positions in their own 
cages. The system is locally rigid and can resist stress. 
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On the much slower time scale Tijje) ^ ^r) 1 ' ^-/^FiHe), 
the cages start to drift around each other leading to a 
plastic deformation expected in a liquid. Here E F ^ He ) 
is the Fermi energy of liquid 3 He. The characteristic 
time TtHe) is associated with a collective rearrangement 
of "equilibrium" particle configurations due to many- 
particle quantum tunneling |^,^|. It determines all 
low temperature Fermi liquid parameters of the sys- 
tems. At t > T( ffe ) the quantum-mechanical indistin- 
guishability of particles plays a decisive role in deter- 
mining the ground state of the system. The energy un- 
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termines the order of magnitude of the temperature of 
the quantum degeneracy E^He)- The quasiparticle ef- 
fective mass fnt He \ can be introduced in the usual way 

Ed(He) — PF(He)l^ m *(He)- ^ can be extracted from 
the temperature dependence of the electron heat capac- 
ity. Accordingly, at T = 0, both the effective mass 
m *He) ~ ? lT {He) n ^ m (He) > an d the magnetic susceptibil- 
ity of the liquid Xh(He) ~ n {He) ^ 2 N T( He) /h » Xo(He) are 
significantly enhanced as compared to their values Xo(He) 
and w,(H e ) m a noninteracting liquid |l7|]lS| l. Here /ujv 
and ri(He) are the nuclear Bohr magneton and the con- 
centration respectively. The latter formula corresponds 
to the Curie susceptibility with T = h/Ttjj e \. It reflects 
the fact the particles change their equilibrium positions 
on the time scale t (He) ■ 

According to the picture of the liquid which is nearly 
solid |l5| , not only the linear spin susceptibility, but also 
all nonlinear ones should be large. It means that the com- 
plete polarization of the electron liquid can be achieved at 
relatively low magnetic fields ^nH c ~ Em He) < E F i He y 

Furthermore, near the crystalization point the mag- 
netic susceptibility of the liquid 3 He is much smaller than 
the susceptibility of the solid 3 He. Therefore the energy 
density of the solid decreases with H faster that the en- 
ergy density of the liquid. This means that the magnetic 
field drives the system toward the crystallization Jl5|,|l6| . 
By itself it does not mean that the effective mass of quasi- 
particles rn* He ^ (H) increases with H . However the calcu- 
lations in the framework of the Hubbard model showed 
p6[ that it is the case and m* He ^(H) - m \ He ){^) > 0, 
which means that T(# e ) (H) is an increasing function of 
H. 

Let us now discuss the T-dependence of the viscos- 
ity of the 3 He liquid r]( He )(T). At low temperatures 
T <C Ed(He) the viscosity J7(jj e ) (T) of liquid 3 He is given 
by the Fermi liquid theory rj^ He ^ ~ T~ 2 Q. There 
is, however, a wide temperature interval Ed < T < 
where the liquid is nondegenerate already, but still it is 
strongly correlated. Surprisingly enough the experimen- 
tal data for the T dependence of rj^ He ^ in this tempera- 
ture interval are unavailable. Theoretically this depen- 
dence was considered in |l7 ,|l8| where it was argued that 
in this temperature interval r\ ~ T . 

Let us turn now to the case of the electron liquid. If 



r s » 1 we have £l p > E F , where il p = {Aite 2 n i l 2 /m) 1 / 2 
is the plasma frequency at the wavelength of order of the 
interelectron distance. Thus the discussed above model 
of strongly correlated liquid can be applied to the electron 
liquid as well. Namely, one can view the electron liquid 
as a system which is close to the Wigner crystallization. 
Due to existence of the short range crystalline order effec- 
tive electron mass m* = Tinr and the spin susceptibility 
X = fj, 2 Tn/Ti are enhanced compared to the noninter- 
acting Fermi gas case, where fi is the Bohr magneton. 
The enhancement of both quantities is governed by the 
parameter r 3> f2~ . Furthermore, since zero tempera- 
ture spin susceptibility of the Wigner crystal xs is larger 
than the susceptibility of the Fermi liquid xl, the energy 
density of the Wigner crystal Es(H\\) = Es(0) — XsH 2 
decreases with H\\ faster that the energy density of the 
Fermi liquid E L (H\\) = E L (0) — xlH 2 . This means that 
the magnetic field parallel to the film drives the electron 
system toward the crystallization An estimate 

for the change of the critical electron concentration is 



n c{ H,)-nM= {XS ' XL)Hl >0 



vl - v s 



(1) 



Here vl and vs are chemical potentials of the liquid and 
crystal respectively. By analogy with the theory of 3 He 
we can assume that near the crystalization point the ef- 
fective quasiparticle mass in the liquid m*(H) and the 
parameter t(H) are increasing functions of H. 

Concluding this section we would like to remind that 
the existence of two different energy scales is not a unique 
property of the quantum liquids near the solidification 
point. It has been known since Frenkel J3l| that melting 
of ordinary classical liquids occurs at temperatures much 
smaller than the interaction energy between atoms. As a 
result, following to ifjlfl , the viscosity of classical liquids 
is governed by thermal-activation of the over-the-barriers 
processes and decreases with T exponentially. 

To describe the temperature and the magnetic field 
dependence of the resistance of the system we have to 
consider the electron system in the presence of random 
elastically scattering potential. Below we consider two 
limiting cases: a model of a potential slowly waryng in 
space and a case when the potential is modeled by short- 
range randomly distributed scatterers. 



III. THE CASE OF A SMOOTH SCATTERING 
POTENTIAL. 

Let us consider a model where the fluctuations of the 
external potential of a relatively small amplitude are 
smooth functions of coordinates. If the electron concen- 
tration n is close to the critical n c the system can get split 
into the regions of a Fermi liquid and a Wigner crystal. 
This model can explain the following experimental facts. 

A. The fractions of volume occupied by the Fermi liq- 
uid and the Wigner crystal depend on n and therefore 
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the system should exhibit a percolation type zero tem- 
perature metal-insulator transition as n decreases and 
the area occupied by the Wigner crystal grows. 

B. Since \S ^ XL, the magnetic field parallel to the 
film drives the electron system toward the crystalliza- 
tion [|l5|,0 and the fraction of volume occupied by the 
Wigner crystal increases with increasing Hn . This leads 
to a big positive magnetoresistance as a function of Hn . 
The magnetoresistance should saturate when H\\ > HS 
and the electron Fermi liquid is polarized. The saturation 
field HS ~ Si < Ep decreases with decreasing n. This 

• • M - • rj 

conclusion is in agreement with the experiment [|1 Of] - 

Depending on the value of n and the amplitude of 
the scattering potential at H\\ > HS the system could 
be either in the metallic or in the dielectric regime. If 
n < n c (H^) the system at H\\ > HS is in the insulating 
regime. It exhibits a giant positive magnetoresistance 
which corresponds to decreasing of the localization ra- 
dius with increasing H^ . 

C. Consider the case n < n c (H^) when at H\\ > Hn 
the system is in the insulating phase. The model pre- 
sented above can also explain why in the presence of 
the parallel magnetic field H\\ > Hn the system exhibits 
a big negative magnetoresistance as a function of the 
component of the magnetic field H± perpendicular to 
the film. Indeed at H\\ > HS electron spins are polar- 
ized and the problem of the magnetoresistance is of the 
single particle nature. It was shown that at relatively 
small magnetic field when Ljj ± ^> £, it is dominated by 
an interference of direct tunneling paths p7|-p9[|. Here 
Lh ± = y/ch/eH_i is the magnetic length and £ is the 
localization radius. In this case the magnetoresistance 
is shown to be negative and the H± dependence of the 
resistance corresponds to a correction to the localization 
radius — £(0)) ~ ^(O)!^*. Existence of such an 
interference effect depends on the spin structure of the 
system ground state p8| , ^o| . At H\\ — the tunneling 
is a collective process which involves an interchange of 
positions of electrons with different spins. As a result, 
final states of the system, which correspond to different 
tunneling paths have different spin configurations and, 
therefore, they are orthogonal. As a result, the single 
particle interference mechanism of the negative magne- 
toresistance can be significantly suppressed in the case 

H \\ = °- 

D. The significant increase of the resistance as a func- 
tion of temperature can also be explained naturally in 
the framework of the presented above model as a conse- 
quence of the Pomeranchuk effect. Indeed, if tempera- 
ture is not significantly smaller than the spin exchange 
energy in the Wigner crystal, the spin entropy of the 
Wigner crystal is larger than the entropy of the Fermi 
liquid. This means that the Wigner crystal regions grow 
with increasing temperature. 

E. The Pomeranchuk effect disappears when H\\ > HS 
and electron spins are polarized. In this case entropies of 



both the liquid and the solid are much smaller than the 
spin entropy of the crystal at Hn — 0. This means that in 
the first approximation the areas occupied by the crystal 
and the liquid are T-independent. This explains the fact 
that in the metallic state at H\\ > HS the T-dependence 
of the resistance is much smaller than in the case H\\ = 

F. The question about the origin of the anisotropy of 
the magnetoresistance in the metallic phase is open. The 
fact that the magnetoresistance in the parallel field is 
larger than the magnetoresistance in the perpendicular 
field implies that the orbital effects are insignificant and 
the magnetoresistance is determined by spin magnetiza- 
tion of the system. A possible explanation of the large 
anisotropy of the spin magnetization induced by the mag- 
netic field can be related to the existence of the Rashba 
term in one-particle spectrum of electrons 

ap[cr x n] (2) 

Here n is a unit vector normal to the conducting plane 
and cr is the vector of Pauli matrices and p is the elec- 
tron momentum. A reliable estimate for the value of a 
in Si-MOSFETs is unavailable at the moment. A rough 
estimate [H implies that the magnitude of Eq.2 can be 
just several times less than the electron Fermi energy. We 
would like to note, however, that a) in the presence of 
the electron-electron interaction the degeneracy energy 
is renormalized (Ed < Ep) and b)the electron-electron 
interaction renormalizes the factor a in Eq.2 p0| . It is 
natural to expect that the latter renormalization is deter- 
mined by the parameter rf2 p 3> 1. Both effects lead to in- 
creasing of the relative amplitude of the Rashba term. At 
last, we would like to mention that the spin-orbit inter- 
action can manifest itself also in the in-plane anisotropy 
when the resistance depends on the relative orientation 
of the current and the magnetic field parallel to the plane 

Hi- 

At last we would like to mention that the resistance 
of order h/e 2 is a typical feature of the percolation tran- 
sition. (See for example p8| .) The difference between 
the presented above model and |2S}| is in the many par- 
ticle nature of the transition which, for example, makes 
it sensitive to the parallel magnetic field. 

A difficulty associated with the model of smooth po- 
tential is the following. In the case of a small amplitude 
of the fluctuations of the scattering potential at T = the 
transition should take place at r s = r c s ~ 38, which is sig- 
nificantly larger than the experimental value r c s — 10—20. 
On the other hand, if the amplitude of the potential is 
large, in principle, a percolation type transition can take 
place at any r s . In the latter case, however, the sys- 
tem breaks into three phases: insulator, Wigner crystal 
and the Fermi liquid. One can neglect the difference in 
compressibilities of the Fermi liquid and the Wigner crys- 
tal. As a result, the electrostatic analysis shows that the 
electron density should increase as a square root of the 
distance from the boundary pi) . Then the assumption 
that the liquid-crystal transition takes place at r s = 38 
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leads to a conclusion that the fraction of the area oc- 
cupied by the Wigner crystal is numerically small and, 
therefore, at r s ~ 10 — 20 the discussed above T and H 
dependences of the resistance should not be pronounced. 
A possibility to have r c s less than 38 in the presence of 
short range scattering potential is considered in the next 
section. 

We would also like to mention that even in the frame- 
work of the presented above model the question whether 
there is a well defined transition or just a sharp crossover 
between metallic and insulating phases is oppen. Indeed, 
in the two-dimensional case an arbitrarily small disorder 
destroys the first order phase transition [Q. The ques- 
tion of whether in this case the first order Fermi liquid- 
Wigner crystal transition is transformed into a second 
order one, or it is destroyed completely is open. (See a 
discussion of this question in |j(J.) 

The question about possible localization in the "metal- 
lic phase" in the presence of a weak disorder is also open. 
All existing theories of the localization are perturbative 
in electron-electron interaction. At r s 3> 1 the potential 
energy is much larger than the kinetic energy and these 
theories are not reliable especially in the case r s 3> 1 
when the renormalized mass is big, which means that 
the size of electron packets which carry current should 
be large as well. It may be that the localization length 
in the metallic phase is larger than the sample size. 



IV. SHORT RANGE SCATTERING POTENTIAL. 

In this section we show that a) in the presence of a 
short range potential the Wigner crystallization takes 
place at r s < 38 and b) one can explain the mentioned 
above experiments using the fact that that the electron 
liquid is nearly solid, but without involving the Wigner 
crystal state. 

Suppose that there is an impurity with a short-range 
potential of a radius a ~ n -1 / 2 , which is embedded into 
the metal with r s 3> 1. Let us now discuss the r s and H\\ 
dependences of the two-dimensional cross-section of the 
quasiparticle elastic scattering on impurities A(T = 0) 
at zero temperature. It has been pointed out pTJ that 
due to the existence of the short range crystalline or- 
der in the Fermi liquid with r s >la large characteristic 
length T^l-p-nT 1 ! 2 ^> n -1 / 2 should exists, where the liquid 
behaves like a solid. Since the impurity pins the liquid, 
the value of A should increase with r as well 



A(T = 0) 



TO 

a 3> a. 

TO 



(3) 



Therefore, it is plausible that the "transition" takes place 
when the impurity concentration Ni equals to a critical 
one 



and strongly correlated regions near impurities overlap. 
Thus in disordered samples the "transition" takes place 
at r s <§; 38 and the critical electron concentration de- 
creases with increasing of the amount of disorder. This 
is in agreement with the experimental observation [^]j7]] 
that the critical concentration n c is lower in samples with 
higher mobility. We would like to note that numeri- 
cal simulations in the presence of disorder give a critical 
value of r s for the Fermi liquid- Wigner crystal transition, 
which is significantly smaller than 38 jb|,[l4}. 

As it has been mentioned above both t{H\\ ) and m* 
are increasing functions of Hu . As a result, the zero tem- 
perature resistance of the system 



p(T = 0) 



h NiA(T = 0) 
e 2 r?, 1 / 2 



(5) 



increases with H\\ and saturates at H\\ > Hp 

Temperature corrections to Eq. 1 should be small as 
long as T Td- Since processes of electron-electron scat- 
tering in semiconductors conserve the total momentum, 
the T-dependence of the resistivity of the system can be 
only due to T-dependence of A(T). 

It has been shown p3| , p4|j34| ] that at low temperatures 
the resistance of the system increases linearly with T. 
This effect is due to the temperature dependence of the 
Friedel oscillations induced by impurities. If r s > 1 the 
value of the scattering cross-section turns out to be an 
increasing function of Hu as well |2fi[] . This mechanism, 
in principle, could explain qualitatively T and Hu depen- 
dences of the resistance of the metallic state. It is hard, 
however, to explain in the framework of p^|25| why at 
i?U > Hu the temperature dependence of the resistance 
in the metallic state is much smaller than at H — and 
why the magnetoresistance in the dielectric state with 
respect to H± is negative. It may be, though, that the 
mechanism of the T dependence considered in [p3|-p5 34 
is relevant in the case of GaAs samples. 

Let us now discuss the T-dependence of the resistance 
in the temperature interval Td < T < f2 p . Though in this 
case the liquid is not degenerate, it is strongly correlated. 
Therefore the electron-electron scattering in the liquid is 
very effective and the local equilibrium is reached in a 
short time on a spatial scale of order n -1 / 2 . As a result, 
the flow of the liquid near an impurity can be consid- 
ered in the framework of hydrodynamics. In the two- 
dimensional case the electron liquid exerts a force on an 
impurity given by the Stokes formula F ~ rju/ ln(r? / 'nua) , 
where u and rj are the liquid hydrodynamic velocity and 
viscosity of the electron liquid respectively. In a system 
with a finite concentration of impurities the logarithmic 

factor in the equation for F should be substituted for 

l li 

ln(l/aA i ). Thus the resistance of the system is 



N iV (T) _! 
P( T ) ~ o o In 



1 



N}' 2 a 



(6) 



1 1 TO 2 

lC ~ I 2 ~ 



(4) 



Strictly speaking the hydrodinamic approach is valid 
when a > l ee . Here l ee is the electron-electron mean 
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free path. However, since a dependence in Eq.6 is only 
logarithmic, the estimate Eq.6 is also valid in the case of 
semi-quantum liquid when a ~ l ee ~ n~ 1 / 2 . The loga- 
rithmic factors in Eq.6 is associated with the Stokes para- 
dox in the two-dimensional case |p2] . It is not universal. 
For example, electron-phonon scattering will change it 
significantly. We will neglect this factor estimating the 
resistance of the system by the order of magnitude. 

The essential ingredient of the picture which leads to 
the Stokes formula for F and to Eq.6 is that the hydro- 
dynamic velocity near impurities is significantly reduced 
as compared to its bulk value. This is the reason why 
the conventional description of the two-dimensional elec- 
tron system with the help of the Boltzmann kinetic equa- 
tion containing the electron-electron and the electron- 
impurity scattering integrals |3(| would give a result com- 
pletely different from Eq.6. In the former case the elec- 
tron distribution function is spatially uniform and, there- 
fore, the resistance of the system is independent of l ee and 
proportional to the electron-impurity scattering rate. 

The viscosity of the electron liquid in the semi-classical 
regime can be estimated in a way similar to [|l7 18 . On 
a time scale smaller than r the structure of the low en- 
ergy excitations in the liquid is similar to the structure 
of excitations in glasses At T <C hil the liquid ex- 
citations are two-level systems with the density of states 
per particle vq ~ C -1 . In the case of the electron liq- 
uid U ~ e 2 e~ 1 n 1 / 2 is the typical interelectron interac- 
tion energy. As a result |17t , contrary to the gas case, 
the viscosity of the liquid -q ~ Hnv V 2 /T decreases with 
increasing temperature. Here V is the typical matrix ele- 
ment of the transition between the states in the two level 
systems. Making a natural assumption that V ~ Ml p we 
get 



1 m 



(7) 



(8) 



Thus, we arrive to the conclusion that at T > the 
resistance associated with the impurity scattering 

ft, N« m* Trf 
e z n to 1 

should decrease with increasing T. Therefore, p(T) may 
have a maximum at T ~ T^. At T = and in the "metal" 
regime (NiA 2 (T — 0) <C 1) contributions from different 
impurities into p(T) are independent and 



ft NjA(T = 0) 
p(T = 0) = — m . 



(9) 



Thus at T = Td Eq. 4 matches the zero temperature 
value of the impurity scattering cross-section Eq. 3. 



V. CONCLUSION 

We would like to mention several consequences of the 
model presented above which can be checked experimen- 
tally. 



1. In the metallic phase at r s 3> 1 the thermopower, 
normalized by p(T — 0), should be relatively large due to 
the strong enhancement of the effective mass m* /to 1. 
Also, the thermopower should be a strongly increasing 
function of Hu because of the corresponding increase in 
m*(if||). In principle, large thermopower in a metal can 
also occur due to localized spins present in the sample 
and the Kondo effect associated with them. However, in 
the latter case the thermopower would be a decreasing 
function of H\\ . 

2. The tunneling density of states in the metal- 
lic regime with r s > 1 is reduced by the factor 
m* (Hn) /m 3> 1 compared to that of nonintcracting elec- 
trons. The tunneling density of states will be additionally 
suppressed by the magnetic field parallel to the film. 

3. In the Wigner crystal state, as well as in the liq- 
uid with r s 3> 1, the electron compressibility is negative, 
whereas at r s <C 1 it is positive Q|. The magnetic field 
applied parallel to the plane drives the system toward 
the Wigner crystal state, and thus it should make the 
compressibility of the "metallic state" at 38 S> r s ^> 1 
even more negative. 

4. Following the presented above picture the saturation 
of the magnetoresistance in the metallic phase in the par- 
allel magnetic field takes place at H\\ > H?, when the elec- 
tron Fermi liquid gets polarized. On the other hand the 
Pomeranchuk effect disappears when the Wigner crystal 
is polarized. This means that the significant tempera- 
ture dependence of the metallic state may disappear in a 

much smaller mag netic field ffii ~ - < m. 

11 fi 
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